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¢ About 8% of Lunar Return radiative heating occurs below 9 km/s 
— Based on current models 


¢ Return from lower altitude (e.g. EFT1) is entirely in this speed regime 
* Radiation phenomena not well validated in this speed regime 


Approach 
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onist peesured in EAST Facility 
24” Diameter tubes for low (<0.1 Torr) pressure 
4” Diameter tube for higher (>0.1 Torr) pressure 


¢ Measurement by between 2-4 spectrometers covering 190-1450 nm 


Conditions Measured 
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¢ 51 shots between 7-9 km/s 
— 33 (27 good) on the 24” Tube (0.01, 0.05, 0.14 Torr) 
¢ 15 from 190-500 nm 
¢ 12 from 500-1450 nm 
— 18 (17 good) on the 4” Tube (0.14, 0.30, 0.50, 0.70 Torr) 
¢ All from 190-1450 nm 


¢ Subset of 10 tests selected for further analysis (1 per 
pressure/wavelength/tube diameter combination): 
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¢ Spectra are resolved in wavelength and position behind shock 


Sample Data (500-1450 nm) 
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Shot 38 (Red) - 8.33 km/s, 0.14 torr Shot 38 (IR) - 8.33 km/s, 0.14 torr 
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¢ Spectra are resolved in wavelength and position behind shock 


Non-equilibrium Analysis 
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* (somewhat) arbitrarily assign +2 cm of peak as “non-equilibrium zone” 

¢ Integral of this, divided by tube diameter, is the “non-equilibrium metric” 
¢ Presented as function of wavelength : “spectral non-equilibrium metric” 


Spectral Non-equilibrium Metric 
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Spectral Non-eq Metric 


(W/cm2-sr-um) 


UVviNis Vis/NIR NIR Camera 
a 


190 390 590 790 990 1190 1390 


Wavelength (nm) 


Non-equilibrium metric composite from 4 different spectrometers 


Spectral Non-equilibrium Metric has units of radiance 
— It is equal to the radiance accumulated through the non-equilibrium zone 
if the non-equilibrium region is optically thin 


0.14 Torr Tube-Tube Comparison (190-500 nm) 
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¢ Spectral metric is larger in 4” tube than 24” tube 
¢ Overlap region of spectrometer is consistent 
¢ CN Contamination in 4” Tube 
¢ Velocities differ, optical thickness may differ 
— Check predictions 


DPLR/NEQAIR Comparison (190-500 nm) 


Entry Systems and Technology Division 


= ——— DPLR/NEQAIR, 8.29 km/s - 24" 
—= — DPLR/NEQAIR, 8.33 km/s - 4" 
© 

on 

25 

cox 

BS 

oO 

©. 

(op) 


190 240 290 340 390 440 490 


Wavelength (nm) 


¢ Some increase in radiation predicted at 8.33 km/s 
¢ Increase is sensitive to rate model 
¢ Prediction does not match data 


Ratio of Data 


Tube Disagreement (190-500 nm) Sy 
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— EAST 
— Predicted 


Wavelength (nm) 


Median disagreement : 46% (cf. 16% predicted) 
— Not clear how much of remaining 30% is due to errors in prediction or 
experiment 
Divergence at low wavelength 
— 24” Tube calibration suspect based on S/N 


CN contamination radiance 


0.14 Torr Tube-Tube Comparison (500-890 nm) Sy 
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¢ Molecular emission (500-700 nm) 
— 4” Tube 30% larger than 24” Tube 

¢ Atomic radiation signifcantly higher in 4” Tube 
— Lines may be optically thick 


Predicted Non-equilibrium metric 
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Spectral Non-eg Metric 
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DPLR/NEQAIR prediction shows larger metric in 4” Tube 
— Indicates atomic lines are optically thick 
Molecular radiation not predicted by NEQAIR 
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Ratio observed in EAST matches predicted ratio for atoms 


Predictive Modeling Sy 
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DPLR/NEQAIR are used to produce 1D (stag. line) profiles for 
comparison to shock tube data 
Three “heritage” modeling options discussed 
— Park90 with Te=Tt (DPLR Default) 
— Park93 with Te=Tv 
— Johnston14 with Te=Tv (LAURA default) 
Revisions to Model will be discussed 
Use data to guide reasonable modeling assumptions 
Use third party measurements of input parameters 
Do not “tune to fit” 
Maintains some level of independence between model and data set 


Spectral Non-eq Metric 
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Analysis will be divided by spectral features for discussion 
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NO Radiance 
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¢ NO Radiance from (primarily) y, « bands 
— Originate from A2x and D2z states 


¢ Also 6 band (C2IT) 


NO Comparison to Heritage 


Entry Systems and Technology Division 


8.18 km/s, 0.01 Torr 7.34 km/s, 0.70 Torr 

1.0 

0.9 —— EAST — EAST 
g 08 — Park93 é — Park93 
2 —~ 0.7 — 3 
ze 06 — sua = paca cine 
eo = 8 — Park90 (Te=Tt) 
ow 0.5 = 
Zc Ss} 
wel 04 = 
= © 
5 = 03 S 
Q oO 
D 02 D 

0.1 

0.0 

190 210 230 250 270 
Wavelength (nm) Wavelength (nm) 


¢ Underpredicted at all conditions, by all models 


NO Boltzmann 
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8.18 km/s, 0.01 Torr 7.34 km/s, 0.70 Torr 
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¢ Boltzmann Radiance is typically an upper bound for non-equilibrium 
radiation (in compression 
¢ Park models cannot match Boltzmann radiance at 0.7 Torr 
— Must check reaction rates 
¢ Boltzmann radiation too high at 0.01 Torr 
— Non-Boltzmann model needs examination 


NO Reaction Kinetics 
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Rate Coefficient (cm~ /mol-s) 
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NO Formation is driven by so-called Zel’dovich exchange Reactions: 
N,+O<~NO+N 
0,+N<@NO+O 
NO Destruction depends on direct dissociation: 
NO+M—~N+0O-4+M 
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We opt to carry rates from combustion literature (Tsang/Baulch) 
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Impact on NO concentration (0.7 Torr) 
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Updating Exchange Reactions increases peak NO density 
Reducing dissociation rate reduces decay 
Changing the ratio of dissociation by atoms vs. molecules further 
increases NO density 

— Johnston follows Park : ratio is 22 

— Figure shows ratio of 1.0 

— Tsang recommended ratio of <1 


NO Non-Boltzmann modeling Sy 
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For these conditions, NO non-Boltzmann is dominated by heavy 
particle processes 


Internal excitation: 
NO(X) + M ~— NO(A,C,D) +M 
Heavy particle impact Dissociation: 
NO(A,C,D)+MoN+O+M 


Internal excitation rates in NEQAIR are only approximate, 
fundamental data is not available 


The reverse of internal excitation is quenching : rates are available 
at 300K. Assume: 
T(K) 
k, =k,¢.|—— 
7 #300 
Heavy particle impact dissociation is updated to be consistent with 
rate chemistry 


Ratio of atomic to molecular driven dissociation is still 
undetermined 


Adjust Atom/Molecule Rates 
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Rates adjusted consistently in DPLR and NEQAIR 
Ratio of 5 matches 0.7 Torr data 
Also matches NO y at 0.01 Torr 


NO 6 is overpredicted at 0.01 Torr 
— Possibly experimental error due to lower sensitivity in this region 


Summary of Model Revisions 
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¢ Flowfield model 


Update NO dissociation and exchange rates to be consistent with combustion 
literature 


— Alter ratio of NO dissociation by atoms vs. molecules to 5 
Electron impact dissociation rate from radiation model used for flowfield 
— Associative lonization controlled by T, 
Update selected charge exchange rates 
¢ Non-Boltzmann Radiation Model - Molecules 
— Heavy particle dissociation rate consistent with flowfield dissociation rate 


— Use quenching rates from literature to calculate heavy particle excitation rates for 
molecules 


— Electron impact dissociation calculation corrected 
¢ Estimate and include contributions from excited states 
¢ Non-Boltzmann Radiation Model — Atoms 
— Excitation rates updated to hybrid of Huo (dipole allowed) and Park (unallowed) 
— Include Associative lonization process 


Results — 0.7 Torr, 7.34 km/s (190-500 nm) 
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¢ NO and N.* underpredictions rectified (mostly) 
¢ N, 2"¢ Positive Somewhat Overpredicted 
¢ Reasonable match to temporal trend 
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Results — 0.01 Torr, 8.18 km/s (190-500 nm) Sy 


Entry Systems and Technology Division 


° — EAST 
S 2 Park93 
D a Johnston1 4 
1 wn 
5 a Park90 (Te=Tt) 
re = ——This Work 
Oo — 
oO 
a 
op) 


Wavelength (nm) 


¢ N,* still overpredicted 
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Results — 0.7 Torr, 7.34 km/s (500-890 nm) 
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¢ N, 15t Positive Matched 
¢ Atomic lines nearly matched 
¢ Reasonable match to temporal trend 


Position (cm) 


Spectral Non-eq Metric 
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Results — 0.01 Torr, 8.58 km/s (500-890 nm) 
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Underprediction N, 1S' Positive Matched 
Extra atomic lines eliminated 

Other atomic lines underpredicted 
Temporal trend shows spike at shock front 
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Results — 0.7 Torr, 7.34 km/s (890-1450 nm) 
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¢ Atomic overprediction eliminated, lines that 
are present are reasonably close 


¢ Missing molecular radiation source (TBD) 
¢ Temporal trend looks ok 
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Results — 0.01 Torr, 8.58 km/s (890-1450 nm) 
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¢ Atomic overprediction eliminated 
¢ Integral matches data 


¢ Spike observed at shock front, trend 
otherwise ok 
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Summary Sy 
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¢ Non-equilibrium Radiation Data Measured from 7-9 km/s at 6 
freestream pressures from 0.01-0.70 Torr 


— Comparison across two tubes with different diameter, calibration source 
indicate confidence in data of ~30% (in UV) or better (Vis/NIR) 


— Presentation focuses on highest and lowest pressure ranges 
¢ Agreement to Predictive (DPLR/NEQAIR) Model has been improved 


— Underprediction of N,/NO resolved by changes to rate chemistry, heavy 
particle excitation rates 


— N,* overpredicted at low pressure, revised rate/excitation model fixes 
underprediction at high pressure 


— Predctiion of atomic radiation improved by 
¢ Changing excitation model (high energy states) 
¢ Including associative ionization in non-Boltzmann model (3p states) 


¢ How does your model do? 
https://data.nasa.gov/docs/datasets/aerothermodynamics/EAST/index.html 


Work to go 
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Low pressure overpredictions of 
— N,*: State specific associative ionization? 
— NO, N, : Pre-dissociation rates? 
Missing molecular features in infrared (high pressure) 
Spike in shock front at low pressure 


Underpredicted atomic lines at low pressure 


non-Boltzmann associative ionization model : needs realistic 
statewise rates 
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Backup 
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Identification of features suggests regions for further analysis 
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Reaction Rates 
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¢ There are between up to 23 reactions rates across the 3 models, 11 of 
which have some differences: 


NO+M <> N +04 M increased by Johnston 

N, +O <= NO +N Johnston used rate from Fujita, 2006 

NO+O << O, +N Johnston uses rate from Bose, 1997 

N +0 & NO* +e Updated Park93, Johnston/Park90 same 
= N +N €& No +e Updated Park93, Johnston/Park93 same 
a : a oe 
Nt +N, & NO +N Missing from Park90, Johnston/Park93 
2 same 
rs O.7+0 <& OF +0, Missing from Park90", Johnston/Park93 
v same 
a N, te <& N +N +4 e Differs across all three chemistries 

O, +e & O, +e Missing from Park90/Park93 


* Ao implemented in DPR 


Revised Kinetic Model 
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Controlling 
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data 
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N. Model 


N. Radiance 
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¢ N, Features from 
— 1% Positive System (B°II—A*II) 500-750 nm 
— 2"4 Positive System (C7II—B°I) 280-390 nm 


N. 1S Positive 
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¢ Underpredicted at all conditions 


¢ Bonus Atomic Lines! 


N, 2" Positive 
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¢ Underpredicted at all conditions 


¢ Partly obscured by N,* radiation at 0.01 Torr 


Update to N, QSS 
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¢ Changing NO rates reduced underprediction @ 0.7 Torr 
¢ Introducing N, Quenching rates brought data into overprediction 


¢ Updating electron impact processes obtains near-agreement 
— Slight underprediction of N, 1%t Positive, overprediction of 2"! Positive 


¢ 0.01 Torr data (not shown) now overpredicted in UV, matched in Visible 
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N.* Comparison to Heritage 
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¢ Underpredicted at high pressure 


¢ Overpredicted at low pressure 
— Park90 gets right magnitude, but transient (not shown) is incorrect 


N.* after updates 
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¢ Discrepancy at higher pressure mostly solved by revisions to rate model 
¢ Low pressure discrepancy remains 


Low Pressure N,* : Controlling Reaction Sy 
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¢ N,* primarily formed by associative ionization: 
N+N<—N, +e 


¢ This rate typically controlled by T, : becomes rapid when thermal 
non-equilibrium is siqnificant 


¢ However, ground state N 
does not cross N.* states 


¢ Reactions proceed through 
metastable (and possibly 
excited) N atoms 


¢ This creates dependence 
on T, 
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Change Controlling Temperature 
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¢ Experimental Radiation profile matches N,* density when T, 
controlling 


¢ The predicted radiance (and profile) does not match, however 
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¢ Atomic Radiation 
— 3p states 
— 3d states 
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N, O 3p Comparison to Heritage 
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* Oatom: 

— 777 nm underpredicted at all cases 

— 845 nm underpedicted high pressure, matched low pressure 
* Natom: 

— Low pressure : Fair agreement 


— High pressure : adjusting for baseline, matched by Park93/Johnston, overpredicted 
by Park90 


N, O 3d Comparison to Heritage 
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¢ Significant overprediction, all lines/pressures 


Internal Excitation Rates 
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Park rates place 3d states at Boltzmann level (overpredicted) 
Huo rates equilibrate all states closer to ionization level 


Zatsarinny rates place highest states near ionization limit, lower 
states progress toward Boltzmann 


Hybrid Huo/Park equilibrates between Boltzman/Saha 


Impact of Excitation Rate on Radiance 
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¢ Revised rates underpredict 3p atomic lines 
¢ Three alternatives eliminate 3d overprediction 
¢ Huo/Park slightly higher than Huo or Zatsarinny 
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Peak Radiance 
7.34 km/s, 0.7 Torr 
= 10,598K 


10,645K 
1.27 X10!’ cm? 


= 2.42 X 10!4 cm? 


Traditionally, QSS balances internal excitation with ionization 
But, lonization accounts for 0.15% of N atom chemistry 
N atom mass derivative is: 


81% exchange reactions 
10% molecular dissociation 
9% associative ionization 


Including Dissociative Recombination in QSS 
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¢ State-wise associative ionization rates assumed proportional to overall 
associative ionization rates 

¢ Preference factors dictate which atomic states are formed from a given ion 
state 

- Best agreement uses literature data for ground state preference, no 
preference for other states of N,* 
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Flip-through of Non-equilibrium Metric 
Comparisons 


Non-equilibrium — 190-500 nm (0.01 Torr, 8.2 km/s) Sy 
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¢ All models underpredict NO 
- N,* overpredicted by T,=T, options, Heritage does ok 
¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.05 Torr, 8.6 km/s) Sy 
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¢ NO still underpredicted 
¢ N,* improving for T,=T, options, Heritage now too low 
¢ N, 2"9 Positive still underpredicted 


Non-equilibrium — 190-500 nm (0.14 Torr, 8.3 km/s) Sy 
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¢ NO still underpredicted 
- N,* slightly over for T,=T, options, Heritage underpredicts 
¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.14 Torr, 8.3 km/s) Sy 
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¢ NO underpredicted 

¢ N,* matched for T,=T, options, Heritage underpredicts 
— CNcontamination accounts for disagreement at 388 nm 

¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.30 Torr, 8.1 km/s) Sy 
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¢ NO underpredicted 

¢ N,* matched for T,=T, options, Heritage underpredicts 
— CN contamination accounts for disagreement at 388 nm 

¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.50 Torr, 7.7 km/s) Sy 
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¢ NO still underpredicted 
¢ N.* being underpredicted 
— Worse for Heritage 
¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.70 Torr, 7.3 km/s) Sy 
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¢ NO still underpredicted 


¢ N,* more underpredicted 
— Heritage and newer models becoming more similar 


¢ N, 2"¢ Positive underpredicted 


Non-equilibrium — 190-500 nm (0.70 Torr, 7.3 km/s) Sy 
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10 cm tube — with Boltzmann state populations 
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¢ NO matched with Boltzmann distribution for Johnston rates 
¢ N,* and N, are overpredicted by Boltzmann model 


Summary 190-500 nm Sy 
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NO is always underpredicted 
N2 2"4 Positive always underpredicted 


N2+ 15t Negative underpredicted at high pressure, overpredicted at 
low pressure 


Non-equilibrium — 500-890 nm (0.01 Torr, 8.6 km/s) Sy 
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¢ Broad features due to N, 1S Positive absent from prediction 

¢ High level (4d,5s) N and O lines absent from data 

* O3p: 777 underpredicted, 845 underpredicted 

¢ N 3p: overpredicted 

¢ Errors cancel out when integrated — radiance appears well matched 


Non-equilibrium — 500-890 nm (0.05 Torr, 8.9 km/s) Sy 
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¢ Broad features due to N, 1° Positive still absent 

¢ High level (4d,5s) N and O lines still overpredicted 

¢ O 3p: underpredicted, but closer than before 

¢ N 3p: matched by Park90/Park93, overpredicted Johnston 

¢ Errors cancel out when integrated — Johnston appears to matched 


60 cm tube N, 1‘t Positive O 3p—3s 
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Spectral Non-eq Metric 
(W/cm2-sr-um) 


Non-equilibrium — 500-890 nm (0.14 Torr, 8.4 km/s) Sy 
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¢ Broad features due to N, 15 Positive still absent 

¢ High level (4d,5s) N and O lines still overpredicted 

¢ O 3p: matched by heritage model, underpredicted other models 
¢ N 3p: overpredicted by heritage, matched other models 


Non-equilibrium — 500-890 nm (0.14 Torr, 8.3 km/s) Sy 
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¢ Broad features due to N, 15 Positive still absent 

¢ High level (4d,5s) N and O lines overpredicted 

¢ O 3p: matched by heritage model, underpredicted other models 
¢ N 3p: overpredicted by heritage, matched other models 


Non-equilibrium — 500-890 nm (0.30 Torr, 8.1 km/s) Sy 
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¢ Broad features due to N, 15 Positive still absent 

¢ High level (4d,5s) N and O lines overpredicted, but less significantly 
¢ O 3p: matched by heritage model, underpredicted other models 

¢ N 3p: further overpredicted by heritage, matched other models 


Non-equilibrium — 500-890 nm (0.50 Torr, 7.7 km/s) Sy 
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Spectral Non-eq Metric 
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¢ Broad features due to N, 15 Positive still absent 

¢ High level (4d,5s) N and O lines overpredicted 

¢ O 3p: matched by heritage model, underpredicted other models 
¢ N 3p: overpredicted by heritage, matched other models 


Non-equilibrium — 500-890 nm (0.70 Torr, 7.3 km/s) Sy 
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¢ Broad features due to N, 15 Positive still absent 
¢ High level (4d,5s) N and O lines overpredicted 
¢ O 3p: underpredicted all models 
¢ N 3p: overpredicted by heritage, matched other models 
— Apparent disagreement due to missing underlying N, radiation 


Non-equilibrium — 500-890 nm (0.70 Torr, 7.3 km/s) 
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N, 1° Positive (B>II—A°?TT) 


¢ Boltzmann matches N, 1° Positive (Heritage slightly over) 

¢ High level (4d,5s) N and O lines overpredicted by Boltzmann 
¢ O 3p matched by Boltzmann (all models) 

¢ N 3p: slightly overpredicted at Boltzmann 


Impact of Alternate N Atom Excitation cross-section aap) 
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¢ Huo excitation cross-sections 
— Eliminate spurious radiation from N 4d, 5s 
— Underpredict N 3p features 


Summary 500-890 nm Sy 
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N, is always underpredicted 
Spurious N and O lines originating from 4d, 5s states 
N 3p lines 
— Matched by Park90 (Te=Tt) at 0.05 Torr, overpredicted elsewhere 
— Matched by Johnston at 0.14-0.7 Torr, overpredicted at lower pressure 
— Matched by Park93 at 0.05-0.7 Torr, overpredicted at lower pressure 
O 3p lines 
— Underpredicted by Park93/Johnston, except at 0.01 Torr 
¢ 845 nm line overpredicted at 0.01 Torr 
— Heritage approach 
¢ Nearly matches 845 nm line from 0.01-0.50 Torr 
¢ Underpredicts 777 nm line, but not badly 
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¢ All lines in this range overpredicted 


Non-equilibrium — 890-1450 nm (0.05 Torr, 8.9 mis) a 
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¢ Most lines overpredicted 
— Park90 matches 1362 nm line 
— N 3p line (939 nm) less overpredicted than others 


Non-equilibrium — 890-1450 nm (0.14 Torr, 8.4 mis) a 
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¢ All lines overpredicted 


Non-equilibrium — 890-1450 nm (0.14 Torr, 8.4 cis) a 
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¢ All lines overpredicted 


Non-equilibrium — 890-1450 nm (0.30 Torr, 8.1 mis) a 
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¢ All lines overpredicted 
¢ N 3p line (939 nm) near match by Park93/Johnston 


Non-equilibrium — 890-1450 nm (0.50 Torr, 7.7 mis) a 
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¢ Most lines overpredicted 
¢ N 3p line (939 nm) matched by Park93/Johnston 


Non-equilibrium — 890-1450 nm (0.70 Torr, 7.3 mis) a 
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¢ Most lines overpredicted 
¢ N 3p line (939 nm) matched by Park93/Johnston 
¢ Continuum (N, Band) not predicted 


Alternate N Excitation Cross Sections 
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¢ Alternate cross-sections underpredict N 3p line 
¢ Other lines near noise limit 
¢ O atoms unchanged 


Non-equilibrium — 890-1450 nm (0.70 Torr, 7.3 mis) a 
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¢ Boltzmann improves background agreement, lines still too intense 


Summary 890-1450 nm Sy 
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¢ Atomic Lines originating from higher states generally over predicted 
¢ One N 3p line is matched well by Park/Johnston from 0.3-0.7 Torr 
¢ Molecular radiation at 0.7 Torr mostly matched under Boltzmann 


Predictive Summary Sy 
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¢ Agreement to Predictive (DPLR/NEQAIR) Model is mixed 
Molecular radiation from N,/NO is underpredicted 


¢ Boltzmann distribution takes up underprediction for N. B state and 
NO radiation 


¢ N, C state is overpredicted by Boltzmann 
— N,* radiation prediction varies with pressure 
¢ At low pressure: overpredicted for T,=T,, matched by heritage model 
¢ Reasonably matched for intermediate pressure range 
¢ Underpedicted at high pressure 
— High lying N, O state radiation overpredicted 
— Radiation from 3p states of N predicted well, except at lowest pressure 
— Radiation from 3p states of O mostly underpredicted 
¢ How does your model do? 
https://data.nasa.gov/docs/datasets/aerothermodynamics/EAST/index.html 
(Test 59 - available soon) 


